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Sexual reproduction presents significant challenges to for- 
mal treatment of evolutionary processes. A starting point for 
systematic treatments of ecological and evolutionary phenom- 
ena has been provided by the gene centered view of evolution 
which assigns effective fitness to each allele instead of each 
organism. The gene centered view can be formalized as a 
dynamic mean field approximation applied to genes in repro- 
duction / selection dynamics. We show that the gene cen- 
tered view breaks down for symmetry breaking and pattern 
formation within a population; and show that spatial distri- 
butions of organisms with local mating neighborhoods in the 
presence of disruptive selection give rise to such symmetry 
breaking and pattern formation in the genetic composition 
of local populations. Global dynamics follows conventional 
coarsening of systems with nonconserved order parameters. 
The results have significant implications for the ecology of 
genetic diversity and species formation. 



The dynamics of evolution can be studied by statisti- 
cal models that reflect properties of general models of the 
statistical dynamics of interacting systems Research 
on this topic can affect the conceptual foundations of 
evolutionary biology, and many appUcations in ecology, 
population biology, and conservation biology. Among the 
central problems is understanding the creation, persis- 
tence, and disappearance of genetic diversity. In this 
paper, we describe a model of sexual reproduction which 
illustrates mean field approaches (the gene-centered view 
of evolution) and the relevance of symmetry breaking and 
pattern formation in spatially distributed populations as 
an example of the breakdown of these approximations. 

Pattern formation in genomic space has been of in- 
creasing interest in theoretical studies of sympatric spe- 
ciation [^^^ . These papers advance our understanding of 
the mechanisms of forming two species from one. How- 
ever, they do not address the fundamental and practical 
problems of genetic diversity and spatial inhomogeneity 
within one species — a population whose evolution con- 
tinues to be coupled by sexual reproduction. Moreover, 
and significantly, these papers do not address the implica- 
tion of symmetry breaking and pattern formation for the 
gene centered view as a fundamental framework of evo- 
lutionary theory. In the following, we demonstrate that 
symmetry breaking and pattern formation invalidate the 
gene centered view (whether or not speciation occurs), 
and that they are important for the spatio-temporal be- 
havior of the genetic composition of sexually reproducing 



populations. This has a wide range of implications for 
ecology, conservation biology, and evolutionary theory. 

Before introducing the complications of sexual repro- 
duction, we start with the simplest iterative model of ex- 
ponential growth of asexually reproducing populations: 



N,{t + \) = \,N,{t) 



(1) 



where Ni is the population of type i and is their fitness. 
If the total population is considered to be normalized, the 
relevant dynamics is only of the proportion of each type, 
then we obtain 



P^{t + l) 



A,; 
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where Pi is the proportion of type i. The addition of 
mutations to the model, Ni{t -|- 1) = XijNj{t), gives 
rise to the quasi-species model [|| which has attracted 
significant attention in the physics community. Recent 
research has focused on such questions as determining 
the rate of environmental change which can be followed 
by evolutionary change. 

Sexual reproduction causes offspring to depend on the 
genetic makeup of two parents. This leads to conceptual 
problems (not just mathematical problems) in evolution- 
ary theory because the offspring of an organism may be 
as different from the parent as organisms it is competing 
against. A partial solution to this problem is recognizing 
that it is sufficient for offspring traits to be correlated 
to parental traits for the principles of evolution to apply. 
However, the gene centered view is a simpler perspective 
in which the genes serve as indivisible units that are pre- 
served from generation to generation In effect, differ- 
ent versions of the gene, i.e. alleles, compete rather than 
organisms. This view simplifies the interplay of selection 
and heredity in sexually reproducing organisms. 

We will show, formally, that the gene centered view 
corresponds to a mean field approximation [lO| ]. This 
clarifies the domain of its applicability and the condi- 
tions in which it should not be applied to understanding 
evolutionary processes in real biological systems. We will 
then describe the breakdown of the gene centered view 
in the case of symmetry breaking and pattern formation 
and its implications for the study of ecological systems. 

It is helpful to explain the gene centered view using the 
"rowers analogy" introduced by Dawkins In this 

analogy boats of mixed English- and German-speaking 
rowers are filled from a common rower pool. Boats com- 
pete in heats and it is assumed that a speed advantage 
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exists for boats with more same-language rowers. The 
successful rowers are then returned to the rower pool for 
the next round. Over time, a predominantly and then 
totally same language rower pool will result. Thus, the 
selection of boats serves, in effect, to select rowers who 
therefore may be considered to be competing against each 
other In order to make the competition between 

rowers precise, an effective fitness can be assigned to a 
rower. We will make explicit the rowers model (in the 
context of genes and sexual reproduction) and demon- 
strate the assignment of fitness to rowers (genes). 

The rowers analogy can be directly realized by consid- 
ering genes with selection in favor of a particular com- 
bination of alleles on genes. Specifically, for two genes, 
after selection, when allele Ai appears in one gene, allele 
Bi must appear on the second gene, and when allele A^i 
appears on the first gene allele -B_i must appear on the 
second gene. We can write these high fitness organisms 
with the notation (1, 1) and (—1, —1), and the organisms 
with lower fitness as (1,-1) and (—1,1). For simplic- 
ity, we assume below that the lower fitness organisms 
are non-reproducing. Models which allow them to repro- 
duce, but with lower probabilities than the high fitness 
organisms, give similar results. 

The assumption of placing rowers into the rower pool 
and taking them out at random is equivalent to assum- 
ing that there are no correlations in reproduction (i.e. 
no correlations in mate pairing) and that there is a suf- 
ficiently dense sampling of genomic combinations by the 
population (in this case only a few possibilities). Then 
the offspring genetic makeup can be written as a prod- 
uct of the probability of each allele in the parent popu- 
lation. This assumption describes a "panmictic popula- 
tion" which forms the core of the gene centered view often 
used in population biology. The assumption that the off- 
spring genotype frequencies can be written as a product 
of the parent allele frequencies is a dynamic form of the 
usual mean field approximation neglect of correlations 
in interacting statistical systems [|l^ . While the explicit 
dynamics of this system is not like the usual treatment 
of mean-field theory, e.g. in the Ising model, many of the 
implications are analogous. 

In our case, the reproducing parents (either (1,1) or 
(—1,-1)) must contain the same proportion of the cor- 
related alleles {Ai and Bi) so that p{t) can represent the 
proportion of either Ai or Bi and 1 — p{t) can represent 
the proportion of either A_i or B-i. The reproduction 
equation specifying the offspring (before selection) for the 
gene pool model are: 

Pi,i(t + l)=p(0' (3) 
Pi,_i(t + 1) = P„i,i(t + 1) = p{t){l-p{t)) (4) 

P_i^_i(t + l) = (l-MO)' (5) 

where Pi i is the proportion of (1, 1) among the offspring, 
and similarly for the other cases. 

The proportion of the alleles in generation i -1-1 is given 
by the selected organisms. Since the less fit organisms 



(1, —1) and (—1, 1) do not survive this is given by p{t + 
1) = Pl^^[t + 1) + P[ _^{t+1) = P{^^{t+l), where primes 
indicate the proportion of the selected organisms. Thus 



P{t + 1) = 



Pi,i(t + i) + p_i,_i(t + i) 



This gives the update equation: 



Pit+l) = 



p{tf 



p{tY + {i^p{t)f 



(6) 



(7) 



There are two stable states of the population with all 
organisms (1,1) or all organisms (—1,-1). If we start 
with exactly 50% of each allele, then there is an unsta- 
ble steady state. In every generation 50% of the organ- 
isms reproduce and 50% do not. Any small bias in the 
proportion of one or the other will cause there to be 
progressively more of one type over the other, and the 
population will eventually have only one set of alleles. 
This problem is reminiscent of an Ising ferromagnet at 
low temperature: A statistically biased initial condition 
leads to alignment. 

This model can be reinterpreted by assigning a mean 
fitness (analogous to a mean field) to each allele as in Eq. 
(^. The fitness coefficient for allele Ai or Bi is Ai = p{t) 
with the corresponding A_i = 1 — Ai. The assignment 
of a fitness to an allele reflects the gene centered view. 
The explicit dependence on the population composition 
(an Engligh-speaking rower in a predominantly English- 
speaking rower pool has higher fitness than one in a pre- 
dominantly German-speaking rower pool) has been ob- 
jected to on grounds of biological appropriateness Q. 
For our purposes, we recognize this dependence as the 
natural outcome of a mean field approximation. 

We can describe more specifically the relationship be- 
tween this picture and the mean field approximation by 
recognizing that the assumptions of no correlations in re- 
production, a random mating pattern of parents, is the 
same as a long-range interaction in an Ising model. If 
there is a spatial distribution of organisms with mating 
correlated by spatial location and fluctuations so that the 
starting population has more of the alleles represented by 
1 in one region and more of the alleles represented by —1 
in another region, then patches of organisms that have 
predominantly (1, 1) or (—1, —1) form after several gen- 
erations. This symmetry breaking, like in a ferromag- 
net, is the usual breakdown of the mean field approx- 
imation. Here, it creates correlations / patterns in the 
genetic makeup of the population. When correlations be- 
come significant then the species has two types, though 
they are still able to cross-mate and are doing so at the 
boundaries of the patches. Thus the gene centered view 
breaks down when multiple organism types form. 

Understanding the spatial distribution of organism 
genotype is a central problem in ecology and conserva- 
tion biology |l^,|l^. The spatial patterns that can arise 
from spontaneous symmetry breaking through sexual re- 
production, as implied by the analogy with other models. 
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may be relevant. A systematic study of the relevance of 
symmetry breaking to ecological systems begins from a 
study of spatially distributed versions of the model just 
described. This model is a simplest model of disruptive 
selection, which corresponds to selection in favor of two 
genotypes whose hybrids are less viable. Assuming over- 
lapping local reproduction neighborhoods, called demes, 
the relevant equations are: 



pix,t + l) = Dipix,t)) 
Dip) = 



P 



p2 + (l-p)2 



Nr 



\xj\<R 

Nr=\{xj I \xj\ <R]\ 



(8) 
(9) 

(10) 
(11) 



where the organisms are distributed over a two- 
dimensional grid and the local genotype averaging is per- 
formed over a preselected range of grid cells around the 
central cell. Under these conditions the organisms lo- 
cally tend to assume one or the other type. In contrast 
to conventional insights in ecology and population biol- 
ogy, there is no need for either complete separation of or- 
ganisms or environmental variations to lead to spatially 
varying genotypes. However, because the organisms are 
not physically isolated from each other, the boundaries 
between neighboring domains will move, and the domains 
will follow conventional coarsening behavior for systems 
with non-conserved order parameters. 

A simulation of this model starting from random ini- 
tial conditions is shown in Fig. This initial condition 
can arise when selection becomes disruptive after being 
non-disruptive due to environmental change. The forma- 
tion of domains of the two different types that progres- 
sively coarsen over time can be seen. While the evolu- 
tionary dynamics describing the local process of organ- 
ism selection is different, the spatial dynamics of domains 
is equivalent to the process of coarsening / pattern for- 
mation that occurs in many other systems such as an 
Ising model or similar cellular automata models Jl7| , p^ . 
Fourier transformed power spectra (Figs. ^-^) confirm 
the correspondence to conventional coarsening by show- 
ing that the correlation length grows as t^/'^ after initial 
transients. In a finite sized system, it is possible for one 
type to completely eliminate the other type. However, 
the time scale over which this takes place is much longer 
than the results assuming complete reproductive mixing, 
i.e. the mean field approximation. Since flat boundaries 
do not move except by random perturbations, a non- 
uniform final state is possible. The addition of noise will 
cause slow relaxation of flat boundaries but they can also 
be trapped by quenched (frozen) inhomogeneity. 
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FIG. 1. Spatially distributed evolution with disruptive se- 
lection giving rise to two types appearing in patches and coars- 
ening. The space is periodic and has 256 x 256 sites, and the 
mating neighborhood radius is i? = 5. 




FIG. 2. Fourier power spectra averaged over ten simu- 
lations of evolutionary processes like that shown in Fig. ^ 
(512 X 512 sites and _R = 1). 
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FIG. 3. Temporal behavior of the peak of a Fourier power 
spectrum in the shown case. Top: The peak frequency kp{t) 
which follows approximately t~^^^ . Bottom: The peak power 
S{kp) which follows approximately t^^^. 




FIG. 4. Collapsed version of the Fourier power spectra 
demonstrating the scaling form S{k) = t^^^ f{kt^^^). 

The results have significant implications for ecology of 
genetic diversity and species formation. The observation 
of harlequin distribution patterns of sister forms is gener- 
ally attributed to nonhomogeneities in the environment, 
i.e. that these patterns reflect features of the underly- 



ing habitat (=selective) template. Our results show that 
disruptive selection can give rise to spontaneously self- 
organized patterns of spatial distribution that are inde- 
pendent of underlying habitat structure. At a particular 
time, the history of introduction of disruptive selection 
events would be apparent as a set of overlapping patterns 
of genetic diversity that exist on various spatial scales. 

More specific relevance of these results to the theoreti- 
cal understanding of genetic diversity can be seen in Fig. 
^ where the population averaged time dependence of p 
is shown. The gene centered view / mean field approx- 
imation predicts a rapid homogenization over the entire 
population. The persistence of diversity in simulations 
with symmetry breaking, as compared to its disappear- 
ance in mean field approximation, is significant. Implica- 
tions for experimental tests and methods are also impor- 
tant. Symmetry breaking predicts that when population 
diversity is measured locally, rapid homogenization sim- 
ilar to the mean field prediction will apply, while when 
they are measured over areas significantly larger than the 
expected range of reproduction, extended persistence of 
diversity should be observed. 
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FIG. 5. Comparison of the time dependence of type prob- 
ability in the mean field approximation and symmetry break- 
ing, calculated using different random number sequences. Di- 
versity persists much longer in the latter. In some cases, for- 
ever. 

The divergence of population traits in space studied in 
our work can also couple to processes of speciation, i.e., 
processes that prevent interbreeding or doom the progeny 
of such breedings. These may include assortative mat- 
ing, whereby organism traits inhibit interbreeding. Such 
divergences can potentially lead to the formation of mul- 
tiple species from a single connected species (sympatric 
speciation). By contrast, allopatric speciation, where dis- 
connected populations diverge, has traditionally been the 
more accepted process even though experimental obser- 
vations suggest sympatric speciation is important. 

Recent studies have begun to connect the process 
of symmetry breaking to sympatric speciation. With- 
out considering pattern formation in physical space, we 
and other researchers have been investigating the role of 
pattern formation in genomic space as a mechanism or 
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description of sympatric speciation. These studies in- 
clude: a model of stochastic branching and fixation of 
subpopulations due to genetic drifts and local reproduc- 
tion in genome space general reaction-diffusion Tur- 
ing pattern formation models in genomic space 
and specific individual-based models of reproductive iso- 
lation involving assortative mating and disruptive selec- 
tion (intrinsic disruptive selection, or disruptive selection 
arising from competition or sexual selection) Our 
work, presented here, is unique in discussing spatial in- 
homogeneity and genetic diversity within one species. 

In conclusion, in formalizing sexual reproduction in 
evolutionary theory, we have found fundamental justi- 
fication for rejecting the widespread application of the 
gene centered view. The formal mathematical analysis 
we presented to demonstrate the lack of applicability of 
the gene centered view is an essential step toward devel- 
oping a sound conceptual foundation for evolution. We 
also showed that the gene centered view breaks down 
for species where local mating and disruptive selection 
give rise to symmetry breaking and pattern formation, 
which correspond to genetic inhomogeneity and trait di- 
vergence of subpopulations. The patterns formed un- 
dergo coarsening, following the usual universal spatio- 
temporal scaling behavior. The slow movement of bound- 
aries between types cause long term persistence of genetic 
diversity through the local survival of (partially) incom- 
patible types. This provides a new understanding of the 
development and persistence of spatio-temporal patterns 
of genetic diversity within a single species. 

One should note that the context in which the gene 
centered view breaks down is of profound significance in 
applied aspects of modern ecology and conservation biol- 
ogy. The preservation of endangered species and ecosys- 
tems is currently at risk due to a dramatic decrease in 
their genetic diversity. We have described the implica- 
tions of our results for the experimental observation of 
genetic diversity in endangered species. Our study of spa- 
tial patterns of genetic diversity in populations may also 
help guide the design of conservation areas and human 
directed breeding programs for endangered organisms. 
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